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A clear consensus exists in the German society that renewable energies have to play a dominant role 
in the future German energy supply system. However, many questions are still under discussion; for 
instance the relevance of the different technologies such as photovoltaic systems and wind energy 
converters installed offshore in the North Sea and the Baltic Sea. Also concerns exist about the cost of a 
future energy system mainly based on renewable energies. In order to be able to address the raised issues 
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on a scientifically sound basis we have set up a new simulation model REMod-D (Renewable Energy 
Model-Deutschland) that models the energy balance of the electricity and heat sector including 
all renewable energy converters, storage components and loads for a future German energy system for 
a whole year based on an hourly energy balance. The target energy systems modeled use a high fraction 
up to 100% of renewable energies to cover the electricity and heat demand (heating and hot water). The 
model includes also energy retrofit of buildings as a measure to reduce future heat loads of the building 
sector. A mathematical-numerical optimizer is applied in order to identify system configurations with 
minimal overall annual cost. In this first part of a two-paper series we describe the methodology of the 
REMod-D model and discuss cost and performance values of all included components and in the second 
part we will discuss the results. 

© 2013 Elsevier Ltd. All rights reserved. 
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1. Introduction 


In the summer of 2011 the German parliament approved new 
ambitious targets for the transformation of the German energy 
system. A key strategy towards this transformation mandates a 
planned continuous increase in electricity from renewable energy 
sources: the fraction of total electricity supply from renewable 
energy sources shall increase to 80% in 2050 (35% in 2020; 50% in 
2030; 65% in 2040, respectively).' Second it was decided to reduce 
Germany's primary energy demand to 80%, as compared to 2008 
values, by 2020 and to 50% by 2050 [1]. In addition Germany 
supports the EU objective of reducing greenhouse gas emissions 
by 80-95% by 2050 compared to 1990 [2]. While the long-term 
targets are clearly defined, significant uncertainty still exists on 
which paths to take to achieve these targets as well as on what is 
the optimal mix on the generation, conversion and demand 
side. However, there is almost no doubt that energy savings 
on the demand side, energy efficiency in conversion chains and 
a large fraction of energy generated from renewable sources on the 
production side are the key elements for achieving these long- 
term goals. 

In this paper we describe REMod-D (Renewable Energy Model- 
Deutschland), a newly developed model which performs a sys- 
tematic optimization on two main sectors of the German energy 
system, namely the electricity sector and the heat sector, focusing 
on a future energy system characterized by a high renewable 
energy fraction. In a second paper we present and analyze results 
of the simulations carried out with the model described here [3]. 
Together the electricity and heat sectors cover about 53% of the 
final energy demand and are responsible for about 62% of the 
primary energy consumption (see Section 2). For the model 
development these two sectors were chosen for the following 
reasons: 


- There is an increasingly strong link between the electricity 
sector and the heat sector: on the one hand an increasing 
number of electric heat pump systems are used for heating 


1 Bundesgesetzblatt No. 42, August 4th, 2011, Page 1634 http://www.bundesre 
gierung.de/Content/DE/Artikel/2011/08/2011-08-05-gesetze-energiewende.html. 


buildings and on the other hand an increasing number of 
combined heat and power systems (CHP) in a wide range of 
capacities is employed. 

- Energy-saving building retrofit measures are considered to be a 
key element in increasing energy efficiency and as such it is 
included into the model. 

- The sectors not covered by the model are fuel based energy 
consumption in the mobility sector (private cars, trucks, air 
transportation) and the industry. The long-term future of the 
fossil-fuel based mobility sector is uncertain under today's 
perspective. Various trends are visible and discussed. These 
include battery-powered cars, hydrogen-powered cars using 
fuel cell technology, trucks employing centralized electricity 
using trolley bus technology and many other concepts. A wide 
diversity exists today with regard to using fossil fuels in 
industrial processes. Many processes are available which have 
very different boundary conditions and varying potentials for 
efficiency measures at a variety of costs. Therefore, the mobility 
and industrial processing sectors remain more difficult to 
implement into a coherent energy-economic model based on 
hourly energy balance calculations as used in our model. As of 
now, they are left out of the model discussed here. However, 
their effect on the overall German energy balance is considered 
in Section 2 and briefly discussed again in the second part of 
this paper [3]. 


Various simulation models for national energy systems, which 
consider a high penetration of renewable energies, exist, but only a 
few models are able to combine the supply and demand from both 
the electricity and heat sectors. Many of the energy system models 
developed in the past focus on the electricity sector and its 
development but rather neglect the heat sector or do not cover 
it in detail. For example, the BALMOREL or the SIVAEL model cover 
the electricity sector but only district heating for the heating 
sector [4]. The EMCAS model focuses in detail on the operation of 
power systems [5]. A widely used example that covers both the 
electricity and the heat sector is the EnergyPLAN model developed 
in Denmark which has been applied to a number of different 
scenarios and regions [6-12]. The EnergyPLAN model calculates an 
optimal system configuration regarding heat production from 
solar thermal collectors, industrial CHP, decentralized and 
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centralized CHP units, boilers, heat pumps and thermal storage. 
The electricity production is covered by photovoltaics, on- and 
offshore wind and conventional condensing power plants. Elec- 
tricity storage is implemented by pumped hydro-power storage, 
batteries and the conversion of electricity into hydrogen. The 
optimization is based on hourly energy balances of supply and 
demand and is built to optimize technical regulatory strategies 
such as meeting heat demand or both heat and electricity demand 
[13]. Thus, the EnergyPLAN model helps to visualize the interac- 
tion and cost of heat and electricity producing and consuming 
technologies, but does not offer optimization of total system costs. 
A typical example for a model that focuses on the optimization of 
minimal energy system cost is the MARKAL model developed by 
the International Energy Agency (IEA) [14]. In contrast to the 
EnergyPLAN model, this model is used to optimize system cost and 
to represent the evolution of an energy system over a time period 
of several years. This model was used by the European Commis- 
sion to determine the so called 20-20-20 targets [15]. Other 
models with a comparable aim are, for instance, IKARUS which 
performs optimizations for a time period of up to 40 years in time 
steps of 5 years [16] or PERSUS which focuses on the minimization 
of energy supply system relevant expenditures [17]. These cost 
optimization tools cover the heat sector to a certain extent, but 
there is no detailed description of the related influences like 
building retrofit or alternative future heating technologies for 
centralized and decentralized systems. To our knowledge, there 
is only one tool that focuses on the detailed description of the 
heating and building sector in the context of a complex energy 
system. This model, called Invert [18], determines the influence of 
policy schemes in the field of heating, domestic hot water and 
cooling through a detailed description of this sector. Nevertheless, 
there is no cost optimization for a one year time period with 
hourly time-steps and also the interaction with the electricity 
sector is not treated in detail. We developed the REMod-D model, 
introduced in this paper, to fill the gap between a detailed 
modeling of the building and electricity sector with our complex 
energy system involving both the electricity and the heat sector. 
Further, the total cost optimization of a complex future energy 
system with a high penetration of renewable energy sources is 
carried out in hourly time-steps over a simulation of one year. 


2. Electricity and heat in Germany's national energy balance 


In 2010 the total German primary energy demand amounted to 
3629 TW h (this value and all of the values used in this section are 
based on those provided in [19] unless otherwise indicated). About 
30% (1113 TW h) are attributed to losses and 70% (2516 TW h) are 
used to cover the final energy demand. In Table 2.1 the distribution 
of final energy between the different application sectors - indus- 
try, mobility sector, tertiary sector (commerce, trade and services) 
and private households - and among different areas of application 
is shown. The table also indicates which areas of application are 
covered or not in the REMod-D model. The main areas not covered 
in the model are non-electrically produced process heat in the 
industrial and tertiary sector and non-electrically produced 
mechanical energy in the mobility sector. In both of these areas, 
mainly fossil fuels are used today for either heating of industrial 
processes or in internal combustion engines for traction of private 
cars, trucks and non-electrically driven railed vehicles and for air 
transportation. 

Summarizing for Germany in 2010, the overall heating demand 
for all buildings (private households, tertiary and industrial build- 
ings) was 781 TWh, the hot water demand in all sectors was 
105 TWh and the total electricity demand was 516 TWh. The 
electricity demand was higher in the years before the financial 


Table 2.1 

Distribution of final energy consumption in Germany (2010) for different applica- 
tion sectors and different uses; values highlighted with bold face are covered in the 
REMod-D model. 


2010 in TWh 

Industry Mobility Tertiary Private Total 

households sectors 
Heating 54.4 3.6 187.9 534.9 780.8 
Thereof electricty 0.9 0.8 10 271 38.8 
Hot water 6.3 0 15.9 831 105.3 
Thereof electricty 0 0 3.4 15.6 19.1 
Process heat 463 0 29.4 36.4 528.7 
Thereof electricty 37.5 0 7 35.7 80.2 
Comfort cooling 4.6 0.7 2.8 0 81 
Thereof electricty 4.6 0 2.2 0 6.9 
Industrial 4.9 0 91 25.8 39.8 
refrigeration 
Thereof electricty 4.9 0 9 25.8 39.7 
Mechanical energy 153.5 699.5 59.9 31 916 
Thereof electricty 150.7 14 30.5 31 198.3 
Information 8.8 2.9 21.2 22.3 55.2 
technology 

Thereof electricty 8.8 0.8 21.2 22.3 53.1 
Artificial lighting 10.4 3.5 56.9 11.3 82.1 
Thereof electricty 10.4 0.8 56.9 11.3 79.4 
Total uses 705.9 710.2 382.9 716.9 2516 
Thereof electricty 217.9 16.5 140.1 141 515.5 


crisis. The corresponding values for the total electricity consump- 
tion were 539.6 TW h in 2006, 541.2 TW h in 2007 and 542.2 TWh 
in 2008 respectively [19]. Therefore we decided to base our 
calculations on an electricity demand of 500 TWh without elec- 
tricity used for heating (39TWh) and hot water (19TWh). 
Heating and hot water are directly handled in the model and are 
therefore subtracted from the total electricity demand. The 
remaining electricity loads pertain to all applications except 
heating and hot water. 


3. Modeling approach 


The basic concept of the REMod-D model is to set up a fixed 
topology of energy producers, converters, storage devices and 
consumers and to optimize their sizing to meet the goal of 
minimized annual overall cost. The optimization rescinds - at 
least in this first step - completely from the existing infrastructure 
of the German energy system, i.e., the approach can be character- 
ized as a “Green Field” simulation. 

As a first case, a system has been defined which covers the 
demand in the electricity and heat sector completely by renewable 
energies harvested in Germany and in which no energy exchange 
with neighbor countries is allowed (further called “100% sce- 
nario”). A schematic of the implemented complete system is 
shown in Fig. 3.1. This system contains photovoltaic power 
systems, on- and offshore wind generators and hydro power 
stations as primary electricity generators. Solar thermal collectors, 
connected either to central heating systems or to decentralized 
heating systems, act as primary heat generators. The secondary 
electricity generators are combined heat and power systems 
(centralized, decentralized) and combined cycle power plants. 
Both use either gas from biomass or gas produced from primary 
renewable electricity (power-to-gas). The electric load is given by 
the total load of all electricity consumers connected to the public 
grid which are not used for heating and hot water plus the 
electricity loads that are due to electrically driven heat pumps. 
Heat loads are separately treated for buildings which are con- 
nected to a centralized heat production (district heating networks) 
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Fig. 3.1. Topology of the model of the German energy system (electricity and heat sector) of the REMod-D model. 


and buildings which have a decentralized heating system. 
For buildings having a decentralized heating system, three technolo- 
gies are distinguished, namely electrically driven heat pumps, decen- 
tralized combined heat and power systems and gas driven heat 
pumps. The latter is considered as the follow-up technology of the 
gas and oil boilers (condensing, non-condensing) predominantly used 
today in Germany. Five different types of storage are included in the 
model: pumped hydro-power storage and batteries for the direct 
storage of electricity, gas storage (in Germany typically large caverns) 
for storage of gas from renewables (power-to-gas) as well as centra- 
lized and decentralized heat storage. 

Today two main technologies are being discussed for the 
conversion of renewable electricity into synthetic fuel for long- 
term storage: (1) The production of hydrogen and its use in various 
applications and (2) the conversion of the produced hydrogen into 
methane (Sabatier-process) for use in the existing natural gas 
infrastructure (pipes, storage devices). While the second path has 
the disadvantage of lower overall efficiency and higher converter 
costs due to the additional conversion step, it has the advantage of 
being fully compatible with the existing natural gas infrastructure. 
For our model we therefore chose to investigate the case of 
methane production from renewable energies since the cost for 
all of the involved components that are necessary to maintain a 
natural gas grid are more transparent. In future models we aim to 
integrate the cost of the hydrogen infrastructure and to compare 
the two options. 

A single simulation is carried out as follows: 


© Capacity values are defined for all generation, conversion and 
storage components of the energy system. 

è An hour-by-hour simulation is carried out for a typical year. 
A complete energy balance is computed for each hour of the 
year following the operation scheme described in Section 5. 

e This annual simulation is repeated with modified values of 
the capacity of photovoltaic until the annual energy balance 
is fulfilled, i.e., until long term storage devices (centralized heat 


storage, gas storage) have the same storage content at the 
beginning and end of the year. 

© The total annual cost of this system is calculated. This cost is 
only composed of re-investment in order to replace compo- 
nents at the end of their lifetime and operation and main- 
tenance (O&M) cost. No cost for fossil fuels or uranium for 
nuclear plants arises since the system is operated by 100% 
renewable energy sources (for more details see Section 6). 


Many simulation runs are performed for different parameter 
sets, i.e., different compositions of the capacities of all compo- 
nents. A systematic parameter modification is carried out in order 
to identify the system composition which leads to minimized 
overall annual cost. A fast convergence algorithm has been devel- 
oped and implemented for parameter modification in order to 
achieve a small computer run time and to carry-out sensitivity 
studies in affordable time (for more details see Section 6). 

Beside the “100% scenario”, two further scenarios were con- 
sidered: a scenario in which import of electricity is permitted and 
a scenario which allows the use of some fossil fuel. In these 
scenarios the general approach of simulation and optimization 
is not changed. However, in the import scenario a maximum 
capacity of import of electricity has been fixed and the result on 
the sizing of the system components has been studied. In the 
fossil fuel scenario a pre-defined amount of fossil fuels is used and 
all components are optimized under this modified boundary 
condition. 


4. Models of single components 


In the following, the models used and cost figures for all 
components as well as all assumptions made are described. Cost 
figures are based on expected future costs based on learning curve 
models. This means that for all technologies a future cost value is 
assessed which is expected to be achieved after broad market 
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introduction and exploitation of cost saving potentials by large 
scale production and main steps of technology advancement. 
Today different technologies have, of course, achieved different 
levels of development and thus the distance to the assumed future 
cost value is different. We investigated the important literature on 
each technology and in case of different statements regarding 
expected future cost, we used average values in our model. 


4.1. Load components 


4.1.1. Electric load 

The annual German electricity demand used in the calculations 
amounts to 500 TW h without electricity used for heating and hot 
water. The electricity for heating and hot water is determined by 
the simulation (see Section 4.3 on electrically driven heat pumps). 
This value is assumed to remain constant, i.e., no cost related 
efficiency measures have been implemented in the model and 
thus they have not been subject of optimization. To distribute the 
electricity demand over 8760 h in a year, we used the hourly load 
profile of 2011 published by the European Network of Transmis- 
sion System Operators for Electricity (ENTSOE). Such load data 
is referred to as the hourly average active power absorbed by 
all installations connected to the transmission network or to the 
distribution network. As it is calculated as an average value for 
every hour, it is not entirely precise but of sufficient quality for the 
purpose of this model [20]. 


4.1.2. Heating load 

The annual heating energy for the entire German building stock 
amounts to 780.8 TW h (2010) [19]. This heating demand corre- 
sponds to the total amount of all buildings in the residential, 
tertiary and industrial sector. All these buildings together repre- 
sent a heated useful floor area of 6300 x 10° m? (3400 x 10° m2 
residential [21], 2.300 x 10° m? [22] and ca. 600 x 10° m? based on 
own calculations”). The hourly heating demand, Pheating IS calcu- 
lated in a simple way by the following equation: 


if Tamb > Tset 
if Tamb <Tset 


0 
Pheating = l k(Tset _ Tamb) a) 

Tamb denotes the ambient air temperature (dry bulb tempera- 
ture) of the particular hour and Tet a set temperature for operation 
of heating systems. In the simulation Tset was set as 15°C. k is a 
constant factor chosen such that the overall sum equalizes the 
total annual value. The resulting annual profile of daily heating 
energy sums is displayed in Fig. 4.1. 

In the German energy policy energy-saving building retrofit 
is seen as one of the key measures for reducing the consumption 
of fossil fuels and related CO» emissions (see e.g. [1]). In order 
to include energy-saving building retrofit into the optimization 
process of the REMod-D model, an average cost curve for retrofit 
measures has been developed. A number of studies on the cost of 
energy-saving building retrofit in Germany indicate that the 
specific cost, i.e. investment per unit of saved energy, increase 
disproportionate with an increase in the achieved level of energy 
standard [23-26]. Fig. 4.2 shows the cost curve which was 
approximated based on numbers from [23,25]. It is important 
to note that the cost covered by this function represents only the 
extra value for energy-saving retrofit measures. Basic building 
retrofit to keep a building operational without reducing its energy 
demand is not covered in this cost value. For the purpose of our 
model is assumed that this extra cost is also needed in future 


? This value is estimated based on the total amount of heat for space heating in 
the industrial sector (ca. 54TWh, cf. Table 2.1) assuming that the structure is 
similar to the structure in the tertiary sector. 
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Fig. 4.1. Annual time pattern of daily low temperature heating demand. 
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Fig. 4.2. Average investment cost for energy-saving building retrofit per unit of 
floor area in Germany as used in the simulation. The x-axis denotes the level 
of heating energy demand per floor area in a percentage of the average value of 
2010 (average value 124 kW h/(m?2a) for all buildings). 


Table 4.1 
Summary of used values for electricity, heating and hot water load. 
Load type Annual demand Cost 


Subject of Data sources, 


(2010) in TWh information optimization references 
Electricity 500 n.a. No [19,20] 
Heating 780.8 See Fig. 4.1 Yes [19,21-23,25] 
Hot water 105.3 n.a. No [19] 


building retrofit after the life cycle of energy saving measures has 
been achieved. 


4.1.3. Hot water load 

The annual energy demand for hot water in the total German 
energy sector amounts to 105.3 TW h (2010) [19]. This hot water 
demand corresponds to the total amount of all buildings in the 
residential, tertiary and industrial sector. For simplicity reasons 
this value was equally distributed over all hours of the year, i.e., 
the amount in each single hour is 12,015 MW h. It is assumed that 
this value remains constant, i.e., no cost related efficiency mea- 
sures have been implemented in the model. 

Key values on load components are summarized in Table 4.1. 


4.2. Generation components 


4.2.1. Offshore wind 

To simulate the electricity generation by offshore wind power 
plants, a feed-in curve is used which represents the quarter-hourly 
feed-in of all German offshore plants in the North Sea connected to 
the transmission grid of TSO TenneT in 2011 [27]. This feed-in 
curve has been transformed by building the weighted arithmetic 
mean of every quarter of an hour and normalizing all values on 
the average installed capacity. The resulting curve represents 
roughly 3500 full load hours. Used economical values are 1650 € 
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Fig. 4.3. Electricity generation profiles used in the simulation (daily averages of 
hourly electricity production for the four renewable energy technologies per GW of 
installed capacity). Top: PV and hydro; bottom: offshore wind and onshore wind. 


investment per kW of installed capacity and a lifetime of 25 years 
[28-30]? [31]. The annual profile (daily averages) of electricity 
generation by offshore wind per GW of installed capacity used in 
the simulation is shown in Fig. 4.3. The used values are summar- 
ized in Table 4.2. 


4.2.2. Onshore wind 

For the electricity generation of wind onshore in Germany in 2011, 
the same method as described above is used, applying a different 
source [32]. The feed-in curve has been transformed by normalizing all 
values on the average installed capacity of the corresponding month, 
since continuous installation and re-powering takes place. The used 
values (annual full load hours, specific investment cost and lifetime) 
are summarized in Table 4.2.The annual profile (daily averages) of 
electricity generation by onshore wind per GW of installed capacity 
used in the simulation is shown in Fig. 4.3. 


4.2.3. Photovoltaic systems 

For the electricity generation and the transformation of the feed-in 
curve of photovoltaic in Germany in 2011 the same method as 
described above is used, applying a different source [33]. The used 
values (annual full load hours, specific investment cost and lifetime) 
are summarized in Table 4.2. The annual profile (daily averages) of 
electricity generation by photovoltaic systems per GW of installed 
capacity used in the simulation is shown in Fig. 4.3. 


4.2.4. Hydropower 

To simulate electricity generation by hydropower stations, a 
feed-in curve representing the hourly feed-in of a share of 
installed hydro power stations in Germany in 2011 is used as a 
basis [34]. Due to its incompleteness, the values of this feed-in are 
normalized and then scaled up on the installed capacity. The used 
values (annual full load hours, specific investment and lifetime) 
are summarized in Table 4.2. The annual profile (daily averages) of 


3 All values from [29] are given in USD and are converted into € with an 
exchange rate of 1.3 USD/€. 


electricity generation by hydropower stations per GW of installed 
capacity used in the simulation is shown in Fig. 4.3. 

Annual profiles of the electricity production (daily sums) of 
the four renewable energy production technologies are shown in 
Fig. 4.3. 

Values used for the modeling of electricity generation compo- 
nents are summarized in Table 4.2. 


4.2.5. Solar thermal collectors 

Solar thermal collector systems enable a direct conversion of 
solar energy into heat. These systems may be connected to district 
heating networks or installed in decentralized installations in 
single buildings. Such systems can cover part of the hot water 
load and part of the space heating load. However, due to the 
seasonal mismatch between the space heating load and solar 
gains, the contribution to the heating is relatively small unless 
seasonal storage is used (see Section 4.4). For solar collectors the 
hourly useful collector energy output, Quse,i, is calculated for each 
hour based on the actual hourly average global radiation on the 
solar collector, G;, and the difference between the actual storage 
temperature, Tstor, and the actual hourly average ambient air 
temperature, Tamp, by the following collector equations: 


Quse.i =A; x Gi x Neolli x 3600 s 


Tec Tan iV] F 
Noll, = co Cy x ( Stor = ann) N 
1 


Ncou,i denotes the collector efficiency and co and cı are the 
optical efficiency and the heat loss coefficient of the collector, 
respectively. The+on the squared bracket indicates that the 
expression is only used if the result is greater than zero and is 
otherwise set to zero. The index i represent three different 
orientations (south-east, south, south-west) and two different 
sites which represent southern Germany (Würzburg) and north- 
ern Germany (Braunschweig). The total collector area, Atot, was 
divided into six different sub-areas and for each area the actual 
global radiation on the tilted collector (45° slope) was computed 
for each hour using the software DESIRE [36]. Using the ambient 
air temperature of each location, the collector equations were 
solved and the total useful collector output for each hour, Quse.tot, 
was calculated by summing the six single values for each sub-area, 


with 


6 
Quse.tot = 2 Qusei- (3) 
= 

It was assumed that 50% of the total collector area is located in 
the south and 50% is located in the north of Germany. And it was 
assumed that 30% of the total collector area is oriented south-east, 
40% is oriented south and 30% is oriented south-west. This 
approach was applied in order to take spatial and angular 
distributions of solar thermal collector systems into consideration. 

The values for optical efficiency, co, and heat loss coefficient, c4, 
as well as cost values used in the simulations are summarized in 
Table 4.3. 


4.3. Conversion components 


Following the topology chart in Fig. 3.1, there are five energy 
converters in the modeled energy system. These converters are 
described below and imply the conversion of power-to-gas, gas- 
to-power, gas-to-power-and-heat, power-to-heat, and gas-to-heat. 
Table 4.4 summarizes all performance and cost values used for 
conversion components. 


4.3.1. Power-to-gas 
The production of synthetic gas (such as hydrogen or methane) 
is one of the options to generate chemical fuels from renewable 
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Table 4.2 
Summary of used values for electricity generation components. 


Generation 
technology 
Offshore wind 3500 1650 25 
Onshore wind 1800 1150 25 
Photovoltaic 975 950° 25 
Hydro-power 4500 2400 50 


Full load hours p.a. Specific cost in € per kW Lifetime in years O&M cost in % of investment Subject of 


Data sources, 


optimization references 
3.0 Yes [27,30] 
2.0 Yes [30,32] 
1.0 Yes [30,33] 
2.0 No (fix value 4,7 GW) [28,29,34,35] 


* This number is calculated assuming that 75% of the photovoltaic facilities are installed on rooftops and 25% in utility scale on other areas. 


Table 4.3 
Summary of used values for solar thermal collectors. 


Solar thermal collectors Efficiency values 


€ per m? 
Large fields Optical efficiency co: 0.8 140° 
Decentralized systems Heat loss coefficient cı: 3.0 W/(m? K) 270° 


Specific cost in 


Lifetime O&M cost in % Subject of Data sources, 
in years of investment optimization references 
20 1.4 Yes [28,29] 

20 1.3 Yes [29,38] 


* The value gathered from [29] is given per kW. To calculate the costs in €/m? a conversion factor of 0.7 kW/m? from [37] is used. 


Table 4.4 


Performance and cost values of conversion technologies used in the REMod-D model. 


Conversion technology Conversion efficiency, COP Lifet-ime 
ina 
Value Unit 
Power-to-Gas 0.6 kWeas/kKWe1 17 
Combined cycle 0.65 kWei/kWeas 30 
CHP, centralized 0.65" kWei/kWeas 24 
0.55" kWei/kWeas 
0.35° kWen/kWeas 
CHP, decentralized 0.33 kWe/kKWeas 25 
0.52 kWen/kWeas 
Electric heat pump See Fig. 4.4 kWin/kWe1 20 
Gas heat pump See Fig. 4.4 kKWin/kWoas 20 
Condensing boiler, centralized 0.98 kWen/KWeas 20 
Condensing boiler, decentralized 0.96 kWen/kWeas 20 


O&M cost in % of Specific cost Subject of Data sources, 

invest-ment optimi-zation references 
Value Unit 

3.0 1500 €/kWe1 Yes 40] 

3.0 640 €/kWe1 Yes* 41-43] 

3.0 650° E/kWaı Yes 29,49,50] 

3.0 1400° €/kWer Yes 29,50] 

15 1050 €/kWin Yes 28,51] 

15 800 €/kWin Yes Own assumptions 

1.0 90 €/kWen No 52,53] 

1.0 670 €/kWin No 52,54] 


* The capacity of the combined cycle is determined in such way that in the worst-case hour the required electricity demand is satisfied. 
> values used for electricity only production, no heat use. Today's highest value is 62% [49] and we assume a further increase until 2050. 
© These values are calculated with today's cost from [50] and are reduced with a cost degression factor of 0.6 and 0.75 for decentralized small scale and centralized large 


scale CHPs, respectively, according to the numbers from [29]. 
4 Combined mode: Electricity production and heat use. 


electricity. Additionally for methane it opens the possibility to use 
distribution and storage systems that are available already today. 
In particular in Germany, a lively discussion is ongoing about the 
production of synthetic fuels as an option to overcome seasonal 
mismatches between renewable energy availability and loads. 
Also the pros and cons of using hydrogen directly or its further 
conversion into methane is actively discussed [39]. In our model 
we implemented a general power-to-gas component in a very 
simple way by assuming a constant conversion efficiency of 0.60 
from power to gas. This is the average conversion value for 
electricity into methane, which is highly dependent on the CO, 
source that is needed for the methane production [40]. The 
potential long-term future cost of this technology has a broad 
range of values. Therefore, we took an average value of 1500 € per 
kWelectric at a lifetime of 20 years [40]. 


4.3.2. Combined cycle power plants (gas-to-power) 

Combined cycle power plants are large scale high-efficiency 
units for the conversion of gas (natural gas or synthetic gas from 
renewables) into electricity. Since they are installed with large 
capacity, it is often difficult to make use of the waste heat. Also, 


the amount of waste heat is comparatively small due to the high 
conversion efficiency for electricity production. We used the 
following values: conversion efficiency of 0.65 and long-term cost 
of 640 € per kWetectric at a lifetime of 30 years [41-43]. 


4.3.3. Combined heat and power plants (gas-to-power-and-heat) 
Combined heat and power plants (CHP) are used for a simulta- 
neous production of electricity and heat from gas (natural gas or 
synthetic gas from renewables). In our model we considered two 
different types of CHP sizes, namely large units which are operated 
in combination with district heating networks and small units 
which are installed in a single building. For large systems we 
assumed that they can be either operated in combined mode 
where they produce electricity and heat or in a mode where only 
electricity is produced. The first operation mode has a lower 
electric performance and is only employed if heat can be used 
either by directly feeding the heat load or charging centralized 
heat storage. Large and small systems differ by both conversion 
efficiency and cost. All values used in our model are summarized 
in Table 4.4. Conversion efficiency values for heat, Near, and 
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electric heat pump 


Eehp» COP ohp 


0.0 + r r r r r r r 1 
20% 30% 40% 50% 60% 70% 80% 90% 100% 


heating demand [100 % = 2010 value] 


Fig. 4.4. Dependence of heat pump performance (eenp for electric heat pumps, 
COP;np for gas heat pumps) from building energy standard, expressed as heating 
demand in a percentage of the 2010-value. 


electricity, Electricity, are defined as: 


Electric (4) 


NHeat = Ec and "electricity = Ec 
as as 


with Egas denoting the energy content of the gas and Q and Eetectric 
representing the useful heat output or electric energy output, 
respectively. 


4.3.4. Electrically driven heat pumps (power-to-heat) 

Electrically driven heat pump systems can be used for space 
heating application as well as for production of hot water. For this 
purpose they may use different heat sources such as the ground, 
ground water or ambient air. These systems enable a much more 
efficient conversion of electricity into heat compared to simple heaters 
based on electrical resistances. Therefore, in our model, only heat 
pumps for the conversion of electricity into heat were considered as 
electrically operated heating systems. However, heat rods can be 
utilized for converting surplus electricity to heat (cf. Fig. 5.3). Based 
on long-term monitoring of many heat pump systems throughout 
Germany, detailed figures about their long-term performance exist 
[44,45]. Based on these results, we assumed a conversion factor of 
electricity into heat, Eenp, defined as 


Q 
ia ~ Electric ©) 

We assumed that this factor depends on the energy standard of 
the building sector. This is based on the fact that low temperature 
heating systems can be employed in buildings with a lower 
heating demand. The lower the supply temperature of the heating 
system, the higher is the energy performance of the heat pump. 
The used function is shown in Fig. 4.3; it takes into consideration 
that with decreasing heating demand for buildings the relative 
share of hot water demand increases. The hot water production 
always requires a temperature of about 60 °C and is not dependent 
on the energy standard of the building. Further values used in our 
model (cost, lifetime) are summarized in Table 4.4. 


4.3.5. Gas-driven heat pumps (gas-to-heat) 

The dominating technologies used today to convert fossil fuels 
into heat are oil or gas boilers. At a maximum, these systems 
convert the energy in the fuel into the same amount of low 
temperature heat. An increasing number of systems installed in 
the heating markets are condensing boilers which also make use 
of the condensing heat of water vapor in the exhaust gas of the 
burner. Recently, a number of developments focus on the sub- 
stitution of this heating technology by more efficient systems, i.e. 
systems which better exploit the exergy content of the fuel. First 
field test results indicate coefficient of performance values, COPgnp, 
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Fig. 4.5. Example of the annual time pattern of gas storage filling status. 
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Fig. 4.6. Example of the annual time pattern of the temperature of centralized heat 
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Fig. 5.1. Flowchart for a single simulation with fixed values of all components 
(except capacity of photovoltaic which is adjusted in order to fulfill the condition 
of a neutral energy balance). ATstor denotes the absolute value of the difference 
between the storage temperature of the previous and the recent iteration step. 
Afgas denotes the absolute value of the difference of the filling status of the gas 
storage of the previous and the recent iteration step. 


defined as 


Q 


COP ghp = Eon (6) 
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Determine total heating load P heating (eqn. 1) 


Determine total heat load Pheat = Pheating + Phot water 


Calculate heat load of all buildings using electrical heat pumps: P heatei up = fp * Pheat 


A |Calculate gains of all solar thermal collectors installed in buildings with electrical heat pumps (eqn. 2, eqn. 3) 
5 
= A Case 1: solar gains exceed actual heat load 
= [33 
3 Sb 3 Use excess heat to load decentralized heat storages 
5 2 Storage full | Remaining solar thermal gains are lost 
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5 Case 2: solar gains lower than actual heat load 
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3 A 3 Cover remaining heat loads by storage discharge 
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< 2 Determine amount of heat needed from electrical heat pumps 
3 storage 
empty a . 
Calculate electricity needed to operate electrical heat pumps 


Calculate heat load of all buildings using gas heat pumps: P heat,cup = foup * Pheat 


Calculate gains of all solar thermal collectors installed in buildings with gas heat pumps (eqn. 2, eqn. 3) 


PA 
a 
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a So 3 Use excess heat to load decentralized heat storages 
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i = Storage full Remaining solar thermal gains are lost 
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218 g |Cover remaining heat loads by storage discharge 
= g = 
= 3 

ec Determine amount of heat needed from gas heat pumps 

3 storage 

empt 
py Calculate gas needed to operate gas heat pumps 


Calculate heat load of all buildings using decentralized CHP: P heat acup = facup * Pheat 


Calculate gains of all solar thermal collectors installed in buildings with decentralized CHP (eqn. 2, eqn. 3) 


= 
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Calculate heat load of all buildings using centralized CHP: Pheatccup = fecup * Pheat 


Calculate gains of all solar thermal collectors installed with centralized CHP (eqn. 2, eqn. 3) 
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Fig. 5.2. Approach (operation cascade) for treating heat sub-systems (solar gains =hourly output of the corresponding solar thermal collector). 


of up to 1.4 ([46-48]. As this technology is still in its early infancy, 
a significant increase in performance can be expected and we 
assumed that COPghp-values as shown in Fig. 4.4 will be reached in 


the future. The values still lie significantly below the theoretical 
limit defined by two combined Carnot-cycles operating in opposite 
directions and exploiting the high temperature of the exhaust gas 
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Calculate total electric load Pioag,el,tot (eqn. 13) 


Calculate total electric generation Pyen ctor (eqn. 14) 


Case 1: Electricity generation exceeds electricity loads 


Charge battery storage 
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Charge pumped hydro-power storage (PHPS) 
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sub-system with electrical heat pumps 


sink temperature of heat pumps 


Heat storage has reached maximum heat] 
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Case 2: Electricity generation lower than electricity loads 


Discharge battery storage 


Battery mily Discharge pumped hydro-power storage (PHPS) 
discharged 


PHPS fully 
discharged 


Operate centralized CHP in heat + electricity mode, if heat storage of 
centralized CHP sub-system not fully charged 


electricity load 
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system fully charged 


Operate decentralized CHP, if heat storage of decentralized 
CHP sub-system not fully charged 


electricity generation < 


Heat storage of decentralized storage 
fully charged 


Operate centralized CHP in electricity-only mode 


sufficient to cover electric load 


Electric capacity of centralized CHP not 


Operate combined cycle power plant“ 


Fig. 5.3. Approach (operation cascade) for treating the electricity sub-system for the “100% scenario” ( the capacity of the combined cycle is sized such that in each hour, 


i.e. also in the worst-case hour, the required electricity demand is satisfied. 


of the gas flame. All values used in our model are summarized in 
Table 4.4. 

Table 4.4 summarizes all performance and cost related values 
of the conversion technologies described above and used in the 
REMod-D model. It should be noted that the cost for all heating 
technologies will be reduced later in the calculation with the cost 
of a standard centralized or decentralized heating technology, 
namely a condensing boiler using gas for district heat or for on-site 
heating, respectively. Since these technologies would have to be 
installed anyway, the basic cost of standard reference technologies 
is not considered. However, we assume that in a future energy 
system using predominantly renewable energies, only technolo- 
gies which enable a reasonable thermodynamic efficiency will be 
justified. The given costs include the cost for the heat source and 
its exploitation. 


4.4. Storage components 


Storage components will play a central role in a future energy 
system with a high share of non-dispatchable renewable energy 
sources. Three principle types of storage are contained in the 
REMod-D model: electricity storage (pumped hydro-power sto- 
rage, batteries), gas storage (caverns) and heat storage (centra- 
lized, decentralized). 


4.4.1. Pumped hydro-power storage 

The total installed storage capacity of pumped hydro-power 
storage in Germany is about 40GWh with a maximum total 
power capacity of about 6.7 GW [55]. A number of expansions 
of existing installations or new storage systems are under 


development or in the planning phase. We assumed that a total 
storage capacity of 60 GW h at a power capacity of 10 GW will be 
available in future. This number is not a subject of optimization 
and has been set as a constant in all simulation runs. The overall 
storage conversion efficiency is assumed as 0.8 [56]. All values 
including lifetime and cost values are summarized in Table 4.5. 


4.4.2. Batteries 

Batteries store electricity in an electro-chemical process with 
high efficiency. Both the need for storage of non-dispatchable 
renewable electric energy as well as the interest in electro- 
mobility concepts for private cars stimulated a lot of R&D activities 
worldwide and lead to a strong technology improvement in this 
sector. A number of studies discusses potential future cost and 
performance of battery systems [42,57]. In the REMod-D model we 
used a very simple battery model which describes the perfor- 
mance by a simple single value for the overall storage efficiency. 
Used values including lifetime and cost values are summarized in 
Table 4.5. 


4.4.3. Gas storage 

Germany's capacity for the storage of natural gas is in the range 
of 35 billion m3, which corresponds to about 350 TW h. In our 
model we assumed that this storage capacity is available for 
storage of synthetic gas produced by renewable electricity in 
power-to-gas converters. During simulation we only controlled 
that the needed gas storage capacity which results from optimiza- 
tion does not exceed this limit. Furthermore we assumed that 1% 
of the energy is lost or used for operation [58] of the gas network 
and that there is no extra investment due to the already existing 
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Table 4.5 
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Performance and cost values of storage systems used in the REMod-D model (PHPS: pumped hydro-power storage) 


Storage type Energy performance Cost parameters Subject of Data sources, 
optimization references 
Performance Value Unit Specific Unit Lifetime O&M cost in % of 
parameter cost ina invest-ment 
PHPS Overall efficiency 0.8 kW haı/ 1600 €/kWel 60 1 No [29,56] 
(out/in) KW hei 
Battery 0.9 KW he/ 300 €/kW he 15° 1 Yes [42,57] 
kW he 
Gas 0.99 kW heas/ - - - - No [58] 
kW heas 
Heat (water), Charge efficiency 0.75 kWhiy/ 20 €/m? 40 1 Yes [59] 
centralized kW ħin 
Time constant 180 Days 
Min. temperature 55 T 
Max. temperature 95 g 
Heat (water), Charge efficiency 0.9 KW hin/ Cost contained in 20 1 No [52] 
decentralized kW hin system cost for 
Time constant 72 h heating system 
Min. temperature 55 C 
Max. temperature 95 °C 


* This is equivalent to around 3000 load cycles, assuming that the battery is used 200 days per year. 


infrastructure. Fig. 4.5 shows an example of the annual filling 
status of the gas storage for one of our simulations. The maximum 
capacity of approximately 80 TWh lies significantly below the 
maximum capacity. The main discharging time is in winter and the 
main charging time in summer. The filling status at the end of the 
year is equal to the filling status at the beginning of the year, 
which confirms that the annual energy balance is fulfilled. 


4.4.4. Centralized heat storage 

Centralized heat storage using water as storage medium is an 
interesting option for long-term storage which can help to over- 
come long-term (up to seasonal) mismatches in the heating sector. 
For instance in Denmark many very large heat storage installations 
have been realized in the past years and many more are recently in 
the planning or construction phase, respectively [59]. All these 
storage devices are connected to district heating networks and in 
many installations large solar thermal collector fields have been 
installed in addition to other heat supply systems such as heat 
pumps or combined heat and power plants. 

The physical model used in the REMod-D simulations is a 
simple single-node model. As such, it underestimates the positive 
impact of stratification which is particularly beneficial for the 
performance of solar thermal collectors. The energy balance of the 
heat storage is given by the following differential equation: 


aT stor 


Cstor ge = N charge x P charge = Paischarge — Pioss (7) 
with 
Cstor 
Pıoss = UAstor(T stor zu; Tamb stor) = Tst x (Tstor = Tamb stor) (8) 
or 


In above equations Cstor denotes the total heat capacity of the 
storage, Tstor the storage temperature, charge the efficiency during 
storage charge (accounting for pipe losses), Pcharge the total heat 
power supplied to the storage, Paischarge the total heat power that is 
removed from the storage, Pjoss the storage losses, UAstor the total 
heat loss coefficient, Tampstor the storage ambient temperature 
(ground temperature in case of a water store in the underground) 
and Tstor the time constant of the storage for self-discharge. This 
differential equation has been solved numerically using the Euler 
method. 

Fig. 4.6 shows the example of the annual time pattern of 
the temperature of centralized heat storage (water storage). 


The storage is mainly employed as long-term storage. All used 
values including lifetime and cost are summarized in Table 4.5. 


4.4.5. Decentralized heat storage 

Decentralized heat storage devices have two main functions 
in the energy system of buildings: on the one hand they can 
overcome mismatches between solar thermal energy supply and 
heat loads and on the other hand they can absorb excess renew- 
able electricity either indirectly through heat pumps or directly by 
converting electricity into heat. The last option should, of course, 
only be used if another thermodynamically more favorable use is 
not possible (see Section 6). 

In our model we assumed that 800 | of hot water buffer storage 
per 10kW of heating capacity are installed for decentralized 
heating systems which interact with the electricity grid (e.g. 
electrical heat pump, decentralized combined heat and power 
system). In the case of gas heat pumps, the storage size is not 
linked to the capacity of the heat pump but to the size of a solar 
thermal collector which may be installed in addition to the gas 
heat pump. Here we assume that 801 of buffer storage are 
installed per m? of collector area. Thus, the capacity of decentra- 
lized heat storage is not a subject of optimization but is strictly 
linked to the capacity of the above mentioned decentralized 
heating systems (which all are subjects of optimization). 

The physical model applied for decentralized heat storage 
is similar to the one applied for centralized heat storage and 
is described above. Values of the parameters used for decentra- 
lized heat storage are summarized in Table 4.5. 


4.5. Other aspects regarding the German energy system 


Other components of the overall energy system are the grids 
(electricity, gas, heat) and other resources such as biomass and 
fossil fuels in cases of less than 100% renewable energy. The import 
and export of electricity has a part in a European transnational 
electricity grid and is an option for energy management. 


4.5.1. Electric grids 

The existing German electricity network has a total length 
of 1.78 Mio kilometers of which 1.16 Mio km belong to the low 
voltage level, about 507,210km to the medium voltage level, 
76,279 km the high voltage level and about 35,708 km to the 
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ultra-high voltage level [60]. While more than 70% of the low and 
medium voltage levels are installed underground, the majority of 
the high and ultra-high voltage levels are installed as overhead 
lines. 

To consider the cost for the needed expansion of the electricity 
grid infrastructure, we calculate average costs to simplify the large 
variety of data. The low voltage grid is mainly influenced by the 
expansion of photovoltaic capacity. To calculate the expenses of a 
grid expansion due to an increase in photovoltaic electricity, we 
multiply the installed capacity with an average value of 140 €/kWaı 
and consider a lifetime of 40 years according to data from [61]. The 
cost to cover the increased grid capacity due to the expansion of 
offshore windmills in the North Sea and the Baltic Sea is con- 
sidered in two steps. First, an average cost of 430 € [62] per kWa 
installed capacity of offshore wind power at an average length of 
20 km and a lifetime of 40 years [63] is implemented to consider 
the interconnection lines that are needed to connect the windmills 
to the land-based grid. Second, due to a higher electricity demand 
in southern Germany, we assume that 40% of the electricity 
produced in the North Sea must be transported through high 
voltage lines to the south. Accordingly, the grid capacity of these 
high voltage lines is sized to deliver this 40% share of the installed 
offshore wind capacity. The cost is assumed with 200 €/kWe, [64] 
at an average length of 400 km and a lifetime of 80 years [64]. 

In 2008 the internal consumption of the electricity production 
sector in Germany, including losses and electricity for pumps in 
power plants, amounted to 72.4 TW h at a total production (includ- 
ing import) of 677.3 TW h [19]. This corresponds to an efficiency of 
89%. Since the majority of the electricity consumption in the systems 
considered here is provided by wind converters and photovoltaic 
systems, the net electricity production of each power generator is 
higher because the need for internal pumps, which are necessary for 
the operation of conventional power plants, will be much lower. Also 
a more decentralized production with a higher spatial coincidence of 
production and consumption will lead to shorter distribution dis- 
tances and thus to reduced losses. Therefore we assumed an overall 
efficiency of 95% for the electricity grid. 


4.5.2. Gas networks 

The total length of the existing German natural gas network is 
about 443,000 km. It operates at three different pressure levels. 
In our model we did not calculate any cost related to the O&M of 
the gas network. We assume that the impact of the sizing of gas 
related components is very small with respect to the overall cost of 
O&M of the gas network. 


4.5.3. Heat networks, district heating networks 

District heating networks enable a common operation of 
various heat producers and consumers that are connected through 
a network which is operated by either water or steam as heat 
transfer fluid. Here the integration of large heat storage devices is 
possible, which assists in overcoming mismatches between heat 
production and heat consumption. 

The major cost of district heating networks is due to the 
connections to the building and the heat transfer units in the 
building, while the far-reaching central network contributes less 
to the overall cost. In our analysis we calculated the overall cost of 
district heating networks as a function of the total maximum 
heating capacity (including heating of hot water) covered by 
centralized heat supply. This value is multiplied with a specific 
average cost value of 400 € per kW [65]. 


4.5.4. Biomass 
The total amount of biomass used in Germany for energy 
applications amounted to 311 TWh in 2010 [19]. About 50% of 


this corresponds to solid biomass (e.g. wood, wood chips), 18% to 
liquids (e.g. biodiesel) and 32% to gaseous biomass (e.g. biogas, 
waste gas). From today's perspective, it is very uncertain how 
much biomass will be available for energy purposes in the future. 
Since almost all biomass products can be transformed into high 
temperature heat, we assume that a major part of future biomass 
used for energy purposes is used for industrial process heat 
generation and - in particular in its liquid or gaseous form - for 
the mobility sector. Due to this uncertainty we assumed that 
50 TW h are available for the electricity and heat sector and that it 
is only used in combination with conversion processes with a high 
thermodynamic efficiency, i.e., not for direct production of low 
temperature heat for space heating or hot water. The amount of 
biomass used has been treated as a constant value in all our 
simulations and is not a subject of optimization. The cost has been 
assumed to be 50 € per MW h ([66]—the average value of biogas 
today is 60 €/MW h). 


4.5.5. Import/export of electricity 

In 2010 Germany exported 59.9 TW h of electricity and imported 
42.2 TW h [19]. It may be an interesting option to increase electricity 
import and export with other countries within the European 
electricity network in the future. This will, in particular, be the case 
when overproduction in one country corresponds with shortages in 
other countries (or vice-versa) due to different weather conditions or 
other factors influencing demand or the production by fluctuating 
renewable energy sources. In general, exchange will be of high 
interest if other countries are able to provide large long-term storage 
capacities (e.g. Austria, Norway). However, it is very difficult to assess 
how much future energy systems can rely on the import and export 
of electricity because the underlying market models and price 
policies are unknown. Therefore we implemented the import and 
export of electricity in a simple way in the REMod-D model by 
defining a maximum capacity of electricity that can be imported. In 
cases where the electricity production from renewable sources and 
centralized combined heat and power systems falls below the actual 
electricity needs electricity is imported up to the pre-defined max- 
imum value. Only in cases where still not enough electricity is 
available to satisfy the load combined cycle power plants are 
employed. A sensitivity analysis allows assessing the impact of this 
operation strategy on the optimal sizing of key components and in 
particular on long-term storage components such as the power-to- 
gas, gas storage and combined cycle subsystem. 


4.5.6. Fossil fuels 

To determine the composition of energy systems for the 
electricity and heat sector of Germany when less than 100% is 
supplied by renewable energy sources, the REMod-D model was 
extended to include a fraction of fossil fuels. We assume that fossil 
fuels are only used in combination with conversion processes with 
a high thermodynamic efficiency, i.e., for instance no classical 
boilers for burning fossil fuels for heating and hot water were 
implemented but rather gas heat pumps. We fix the amount of 
fossil fuels at the beginning of a simulation and then let the 
program determine the optimal composition of all converters. 


5. System simulation, operation and control 


The core of the REMod-D model is the calculation of an hourly 
energy balance of the German energy system according to the 
topologies shown in Fig. 3.1. This calculation is repeated for 8760 h 
(or one year). The calculation for the total year is repeated in a first 
loop until the centralized heat storage has a similar temperature 
at the beginning and the end of annual simulation. A second loop 
is implemented in order to determine the required capacity of 
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photovoltaic’ converters in order to assure that the filling of the gas 
storage has the same value at the beginning and the end of the year, i. 
e. the criteria of an annual neutral energy balance is fulfilled. Then the 
total annual cost related to this composition of the energy system is 
calculated (see Section 6). Fig. 5.1 shows a flowchart of such a 
simulation step for a pre-defined set of all parameters, i.e. sizes of all 
components. 

The models of all single components were described in Section 
4. Here we describe the sequence of calculating the hourly energy 
balance depending on the actual conditions. The general concept 
of the operation follows the strategy always to prefer the conver- 
sion which assures highest thermodynamic efficiency. The next 
conversion option with next highest level of thermodynamic 
efficiency is only accepted if the conversion in the previous level 
is no longer possible due to limited capacity of converter or 
storage. 


5.1. Heat sub-systems 


The approach for treating heat sub-systems is explained in 
Fig. 5.2. It has to be noted that four different technical systems are 
installed to cover heat loads of the total building sector, namely 
electrical heat pumps, gas heat pumps, decentralized combined 
heat and power (dCHP) and centralized combined heat and power 
(cCHP). Solar thermal collectors can be used in addition but will 
not be able to supply 100% of the heating and hot water demand 
and are therefore treated separately. In a single simulation the 
share of the total building sector supplied with either technology 
is fixed. The sum of the shares of all four technologies has always 
to be 1: 


Janr +S cue +facue +fecup = 1 (9) 


For solar thermal collectors we treated solar collectors installed 
in combination with centralized CHP and district heating networks 
separately from solar collectors installed in buildings with decen- 
tralized heating systems. Thus, in a single simulation the solar 
thermal collector area installed in combination with district 
heating, Acoıı,centralized has a fixed value and the solar collector 
area installed in all decentralized heating systems, Acotdecentralized, has 
a fixed value. The distribution of the solar thermal collector area 
among the three ensembles of the building sector using decentra- 
lized heating systems is the same as for the heating technologies 
themselves, thus: 


Aco l,decentralized = Acollel HP + Acoll, GHP + Acoll,dcHP (10) 
with 

Aco 1, el HP =fa HP X Acoll,decentralized» a 1) 
Acoll, GHP =f cup X Acoll,decentralized and (12) 
Acott, dcp =f acHp X Acoll,decentralized (13) 


In each sub-system two main cases have to be distinguished 
depending on whether the actual solar heat generation, Pgen,sol, 
exceeds or falls below of the actual heat load of the particular sub- 
system. 


5.1.1. Electricity sub-system 
The approach for treating the electricity sub-system is 
explained in Fig. 5.3. 


4 We used the size of photovoltaic systems for equalizing the annual energy 
balance. Reason is that the technical potential for photovoltaics is less limited than 
that for onshore wind and offshore wind. More details about technical potentials 
for renewable energy sources are described in the second part of this paper [3]. 


Here, the first step is to calculate the total electrical hourly load, 
Pioad,el,torr Which is composed of the electric load by non-energy 
applications, Pioaq,e1, and the electricity needed to operate electri- 
cally driven heat pumps, Pei, eı up: 


Pioad,el,tot = Pioad.el +Paı, el HP (14) 


The next step is to calculate the total electricity generation, Pgen, 
tot, from all renewable energy sources, i.e. offshore wind, Pwina, 
offshor, ONShore wind, Pwind,onshore, Photovoltaics, Ppy and hydro 
power, Phydro, and by decentralized and centralized CHP, Petacup 
and Peiccup, respectively, in case they are needed to cover heat 
loads (see heat sub-system) and thus produce electricity. 


Pgen,tot = Pwind,offshore + Pwina, onshore + Ppv + Phyaro + Pet dcHp + Pei,ccHp 
(15) 


Two main cases have to be distinguished depending whether 
the actual generation, Pgen,tor, exceeds the actual load, Pioag,ei,tor OF 
Peen,tot falls short of Pioad,el,tot- 


6. Optimization methodology and computational 
performance 


In the previous section the simulation of a single step, i.e. with 
given values of the size of all components, was described in detail. 
In this section, the optimization methodology and the correspond- 
ing computation performance are presented. 


6.1. Calculation of annual total cost 


After finishing a single simulation, the corresponding overall 
annual cost, Cannualtot, is determined for an energy system with 
given values of the size of all components. This cost is calculated 
using the following equation: 


N 1 
Cannual,tot = 2 [Ceres x Pj x (+ mi) 
i 


ta 


1 
= Cspec,ref,heat,centr x Pıoad,heat, max ,centr (at mretcons 
ref,centr 


|Cspec,ref,heat,decentr x Pioad. heat, max ‚decentr 


1 
x (a: Myef,decentr 
ref,decentr 


[ 1 
F Qheat,spec,rel x Cesbr x Apuild.tot (+ M ref, esbr 
ref,esbr 


+[Ebiomass X Pbiomass] + [Erossiı x Prossil | 


+ Eel import x Peimpor| = [Eetexport x | (16) 


The meaning of the different lines in Eq. (16) is: 


Line 1: Each summand with index ‘i’ in the sum in the first line 
of Eq. (16) corresponds to the annual cost of each component of 
the energy system described Section 4. These components are: 

e Electricity generation: offshore wind converters, onshore wind 
converters, photovoltaic systems, hydro-power stations. 

© Heat generation: solar thermal collectors connected to district 
heating systems, decentralized solar thermal collectors. 

è Conversion: power-to-gas converters, combined cycle power 
plants, centralized combined heat and power systems, decen- 
tralized combined heat and power systems, electrically driven 
heat pumps, gas driven heat pumps. 

© Storage: batteries, pumped hydro-power storage, centralized 
heat storage devices. 

© Grids: electricity grids, district heating networks. 


Cspeci denotes the specific cost of this component, P; the installed 
capacity (power, energy content or area, respectively), t; the lifetime 
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and m; the O&M cost (percentage of investment cost). All values were 

specified in Section 4. 
Line 2 and line 3: In the second and third line of equation 15, 
the cost for reference heating technologies in both the cen- 
tralized and decentralized case, respectively, is subtracted. This 
cost would be necessary also in a completely fossil fuel based 
heating sector using state-of-the-art technology. Thereby only 
the extra cost of high efficient systems (electric heat pumps, 
gas heat pumps, combined heat and power systems) compared 
to standard heating technology (condensing gas boilers) is 
taken into consideration for calculating the total annual cost. 
Line 4: The fourth line refers to energy saving building retrofit 
(abbreviated ‘esbr ’). Qneat,specrel is the average specific heating 
demand of the building stock given as a percentage of the base 
value (2010) and Cesbr the corresponding specific cost of energy 
saving building retrofit (see Fig. 4.1). Apuild,tot is the total heated 
floor area of the building stock (see Section 4.1). 
Line 5: In line 5 the cost of biomass is calculated. Epiomass is the 
annual energy amount of used biomass and Ppiomass the cost per 
unit. Since there are manifold different production technologies 
to produce gaseous or liquid biomass, no single technology and 
their corresponding cost parameters has been considered but 
an average price of biomass is taken. In a similar way the cost of 
fossil fuel is calculated in the case of scenarios allowing the use 
of fossil fuel. 
Line 6: In line 6 the cost and revenues for import and export 
of electricity is calculated based on the amount of imported 
electricity, Eeimport, and related average cost, Pel,import and the 
amount of exported electricity, Eeı,export and the corresponding 
revenue, Pel,export: 


6.2. Optimization methodology 


The optimization process is based on a modified multi- 
dimensional regula falsi approach. Basics of the underlying mathe- 
matics and numerical solution can be found in standard textbooks 
(e.g. [67]). 

In this approach, a minimum and maximum possible value, 
Vmin and Vmax, has to be pre-defined for each parameter which is 
subject of optimization, i.e. a parameter range is defined. In most 
cases we chose zero as the minimum value. The maximum value 
was chosen to be at the reasonable maximum capacity of that 
particular component. For the overall heating energy demand of 
the building sector we chose 20% of the 2010 level as the 
minimum value. Any lower heating energy demand would be 
unrealistic.. In a first iteration loop for each parameter subject 
to optimization, annual simulations are carried out with two 
values, namely with a first lower value vı and with a second 
larger value v2: 


V1 = Vmin+(1/m) x (Vmax— Vmin) and 
V2 = Vmin + ((M — 1)/M) x (Vmax — Vmin) (17) 
We worked with different values of m, e.g. m=3, M=4 and 
M=10 and compared the results. The smaller the value of M, the 
closer the chosen values of the particular parameter are to the mid 
value of the parameter range. Vice versa, the larger the value of M, 
the closer the chosen values of the particular parameter are to the 
boundary of the parameter range. Combination of N values that 
are subject of optimization leads to 2" single simulations per 
iteration loop. In simulating the “100% scenario”, we have 11 free 
parameters that are subjects of optimization and thus 2!'=2048 
single simulations are carried out. For each of the single simula- 
tions the total annual cost is computed. 


After the first iteration loop has been finished, the parameter 
set leading to the minimum of the total annual cost is determined. 
Then, for each parameter subject to optimization, the parameter 
range is reduced. In case the value v; of parameter i was the one 
which led to the minimal total annual cost, the minimum of the 
parameter range, Vmin, remains unchanged and the new maximum 
of the parameter range is reduced according to: 


Vmin ‚new — Vmin ‚old and 
Vmax ‚new = Vmin ‚old +L/(L+ 1) x (Vmax ‚old — Vmin old) (18) 


In case that value v2 of parameter i was the one which led to 
the minimal total annual cost, the minimum of the parameter 
range, Vmin, is increased in a corresponding way and Vmax remains 
unchanged. Again, we carry out the simulations with different 
values of L, e.g. L=3, L=10 and L=100. The larger L is, the lower 
the reduction of the parameter space is from one iteration loop to 
the next. Thus, with a large value of L, more iteration steps are 
required in order to achieve the same reduction of the parameter 
space achieved with a smaller value of L. If, for instance, L is set 
to 3, the parameter range is reduced by a factor of 2/3 for each 
parameter subject to optimization in each iteration. After K 
iteration loops, the parameter range has been reduced to (2/3)*. 
With K=12 iteration loops, the final parameter range equals about 
0.771% of the initial parameter range for each parameter which is 
subject of optimization. Fig. 6.1 shows an example of the overall 
results of the optimization. Each dot represents one specific cost 
value, calculated with one specific set of parameters for a single 
simulation of a whole year. During the optimization process 
several thousands sets of parameters are calculated and compared 
by means of the specific total annual cost value. It is evident how 
the iteration loops converge and that the differences between 
different parameter sets are marginal in the final iteration loop. 


6.3. Computation performance 


We applied the modified multi-dimensional regula falsi method 
because it proved to be robust. On a normal state-of-the-art PC, 
the complete calculation of a 100% scenario with 11 parameters to 
be optimized and 50 iteration loops (programmed in Delphi’) 
takes about 8 h. One of the results identified during carry-out of 
many simulations with changed numerical parameters (M and L of 
the previous section) and modified parameter ranges for compo- 
nents was that many local minima exist, which lead to similar 
values of the goal function, i.e., total annual cost. More details are 
presented and discussed in the second paper [3]. 


7. Summary and outlook 


In this paper we describe a novel modeling approach for 
simulating the electricity and heat sector of a potential future 
German energy system based on 100% renewable energies (or a 
high penetration thereof). The model uses simple models for the 
single components. The interaction between the different compo- 
nents and all related energy flows are simulated in detail on an 
hourly basis for a complete year. An efficient and robust optimiza- 
tion method has been developed which allows many simulations 
to be performed within a realistic running time. As a result, 
comprehensive sensitivity studies can be carried out and a wide 
range of parameter variations can be implemented. With this tool, 
favorable combinations of components implemented in future 
energy systems using a high share of renewable energies can be 
determined. The simulation results are intended to help decision 
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Fig. 6.1. Example of a complete simulation. Each dot refers to the annual cost of a single simulation corresponding to a set of parameters during the optimization procedure. 


makers in policy and the energy economy sector to identify robust 
scenarios and assess related total annual costs of the future energy 
system. In the second part of this paper we describe first results 
obtained from the model. 

The chosen optimizing algorithm is very robust. Nevertheless, 
it cannot be assured that the absolute minimum of the multi- 
dimensional parameter space is tracked. We carried out various 
simulation runs using different numerical parameters (parameters 
spaces for the different components, numerical parameters L and 
M described in Section 6). The results showed very similar total 
annual costs with different compositions of the component size. 
This indicates that many local minima exist. Interpretation of the 
results leads to the conclusion that various system solutions are 
possible which lead to almost the same total annual cost. Thus, 
other than techno-economic criteria may be taken into considera- 
tion for the design of future energy systems with a very high share 
of renewable energy sources. 

Various plans exist to further expand the broadness and the 
detail in the tool. On the one hand, we plan to include the mobility 
and analyze, for instance, the impact of different concepts of 
private mobility (e.g. battery based versus hydrogen based or a 
mix of both) on the overall energy system. Also the industrial 
sector, (i.e. fuel based industrial processes, which are not yet 
covered in the recent REMod-D model) shall be implemented, at 
least on an annual balance. In particular, the building sector will be 
modeled in more detail e.g. by including building thermal mass as 
a means to store energy and thus enabling load shifting and 
storage for fluctuating renewable energy sources. Also the typol- 
ogy of the building sector shall be represented in more detail, e.g. 
distinguishing between different building types and building 
energy standards and the various cost functions for considering 
the extra costs associated with an energy retrofit. Our aim is that 
the model can be applied in more detail for the building sector and 
used to compare different measures in the different sectors of the 
building stock. 

In order to increase the usefulness of the model not only 
in the description of target compositions of the energy system 
but also pathways to get there, we plan to implement a metho- 
dology which allows to model and analyze intermediate steps 
between the today's situation and the target system. 

Finally the model can be applied to other regions or countries. 
This will be more straightforward for countries or regions which 
have a climate similar to Germany. For regions or countries in 


other climatic zones, adjustments have to be made. For instance, 
regions in the earth's Sun Belt have a much lower heat demand 
and higher cooling demand. Also in these locations, other tech- 
nologies such as concentrating solar power CSP may be an 
important option for electricity production and will have to be 
implemented. 
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